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Abstract 

We investigate the effects of Lifshitz dynamical critical exponent « on a family 
of minimal D = 4 + 1 holographic superconducting models, with a particular focus 
on magnetic phenomena. We see that it is possible to have a consistent Ginzburg- 
Landau approach to holographic superconductivity in a Lifshitz background. By 
following this phenomenological approach we are able to compute a wide array of 
physical quantities. We also calculate the Ginzburg-Landau parameter for different 
condensates, and conclude that in systems with higher dynamical critical exponent, 
vortex formation is more strongly unfavored energetically and exhibit a stronger 
Type I behavior. Finally, following the perturbative approach proposed by Maeda, 
Natsuume and Okamura, we calculate the critical magnetic field of our models for 
different values of a. 
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1 Introduction 


The AdS/CFT correspondence [T] is one of the most important developments in theo¬ 
retical physics in recent years. Because of the very nature of the duality, it also provides 
a promising new way of studying gauge theories in the strongly-coupled regime, where 
the usual perturbative methods fail to apply. The gauge/gravity duality has been used 
to gain insight in a wide variety of physical systems, such as the quark-gluon plasma or 
in condensed matter theory. In this last category, the gauge/gravity duality has been 
fruitfully applied in the study of high-Tc superconductivity, where the usual BCS model 
of Cooper-pair creation ceases to be valid due to the strong interactions between the 
system’s components. Because of these non-trivial interactions, holographic methods 
are hoped to shed some light in the understanding of these systems. This area of study, 
known as holographic superconductivity^ is currently an exciting and very active area of 
research, (for some reviews, see, e.g. [U [31IH [3 [S] ). 

In very general terms, the usual models of holographic superconductivity are built 
around a local gauge group symmetry breaking by one of the component fields in the 
gravity side, where the gravitational solution is an asymptotically AdS charged black 
hole. This symmetry breaking in the gravity side signals the beginning of a supercon¬ 
ducting phase in the dual field theory. (See, e.g. [313] )• It has been found, however, 
that in some condensed matter systems phase transitions are governed by Lifshitz-like 
fixed points. These exhibit the particular anisotropic spacetime scaling symmetry 

t —^ X^t , X —y Xx , (1*1) 

where z is the Lifshitz dynamical critical exponent governing the degree of anisotropy. 
This anisotropy breaks Lorentz invariance and the systems are non-relativistic in nature. 
Therefore, in order to study such field theories holographically, the dual gravitational 
description has to be modified. Indeed, it was found in [S] that these Lifshitz-like fixed 
points can be described by the gravitational dual 

ds^ = -I- X! 

which, for z = 1 reduces back to the usual AdSd +2 metric, but for z ^ 1 satisfies the 
anisotropic scaling (HH). A black hole generalization of this metric was found in m- 

The purpose of this communication is to explore various aspects of holographic 
superconductivity with Lifshitz-like fixed points, with a particular focus on magnetic 
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phenomena. We do this by starting from a minimal bulk model and by studying various 
choices of condensates. More concretely, we want to investigate how the dynamical 
critical exponent z affects our system with respect to its behavior in the isotropic z = 1 
case. Most of the existing research on the subject was realized in D = 4. See, for instance 
[TTl [T^ [T51 [13] . In [T5] , the authors do make an interesting treatment of the D = 5 
case, but have their interest put mainly on studying different kinds of superconductors 
(s-wave, p-wave, soliton) and on the computation of condensation and conductivity. 
Regarding the study of magnetic effects in a Lifshitz background, we note in particular 
[M [H]. In the first reference the authors also treat the D = 5 case, but using a 
different condensate as the ones we will propose, and with a focus on the applicability 
of the matching method. 

Although we will initially consider our minimal model in general dimensions, we 
will focus our attention on the D = 5 case. The reason for this choice of dimension is 
that, as noted in [niin], dimensionality plays an important role in the way external 
magnetic fields act in the dual superconducting system. The standard argument is 
that in a 2 + 1 (ZI = 4) dimensional superconductor an external 3+1 dimensional 
magnetic field will always penetrate the material because the energy needed to expel 
the field scales as the volume, while the energy that the system gains from being in a 
superconducting state scales as the area. This results in the system being a Type II 
superconductor. In the case of a 3 + 1 {D = 5) dimensional system such as the one 
we study, both energies scale as the volume and one has therefore a direct competition 
that does not excludes the possibility of obtaining a Type I superconductor. Also, 
while high-Tc samples are typically composed of 2-dimensional Cu 02 layers {cuprate 
superconductors), it is important to examine the effect of thickness when the system is 
probed by external magnetic fields. 

In this respect, in this paper we see that it is possible to have a consistent Ginzburg- 
Landau phenomenological approach to holographic superconductivity |19] in a Lifshitz 
background. We then apply this Ginzburg-Landau approach to compute, among other 
physical quantities, the Ginzburg-Landau parameter of the system, and to see how it 
is affected by the dynamical critical exponent z. We will also study the effect of an 
external magnetic field acting directly on the system, using the approach proposed in 
[20) . In order to have a more complete study of the system’s properties, we managed to 
study a wide array of condensation cases, always within the D = 5 framework, so that 
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the general tendencies in the behavior of physical quantities become more clear. 

The article is organized as follows. In section 2 we describe our minimal model of 
D = 5 holographic superconductivity with Lifshitz dynamical scaling, and present the 
different cases of condensates we will consider. In section 3, we study small fluctua¬ 
tions around the component fields of our model. From these fluctuations we are able 
to compute the penetration and correlation lengths of the superconductor for our dif¬ 
ferent cases. In section 4 we show that the system can be consistently described in a 
phenomenological Ginzburg-Landau approach. We compute the Ginzburg-Landau pa¬ 
rameter K for the different cases in our model. In section 5, we compute the critical 
magnetic field of the superconductor. Finally, in section 6 we summarize the main 
results and discuss some open problems. 


2 Minimal Holographic Superconductor in D = d-\-2 
Lifshitz Background 

2.1 General Setup 

As mentioned in the introduction, the D=d+2 gravitational dual (HH) can be gener¬ 
alized to a black hole solution m 


ds^ = L" -r^^f{r)dt^ + 


dr^ 




where 


f{r) = 1 - 






( 2 . 1 ) 


( 2 . 2 ) 


and where Vh is the horizon of the black hole. The Lifshitz dynamical critical exponent 
can take values 1 < z < d. The gravitational solution (EH) can be obtained from the 
action m 

® = mikr, / 5 


with extremal solution for the fields 




3^rt — 


\<p _ \y/2{z-l)d 


= 


e ^ = r 

..2 _ r, r 2 


2d 

z — 1 ’ 

(z -I- d - l)(z -I- d) 


q = 2L^{z — l){z + d), 


A= - 


2L2 


(2.4) 
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For the remaining of the paper we will set = 1. We will also prefer to work with the 
coordinate u = Vh/r. This change of coordinates gives 


where 


J2z 


ds^ = — - 


/(w) 


,2z 


dt^ 




du‘‘ 


r: 


d 

IE 

2=1 


dxf 


f{u) = l-u^+\ 


(2.5) 


( 2 . 6 ) 


and the Hawking temperature is 


Th 




(2.7) 


It is therefore the actions (12.3|) that will provide us the gravitational Lifshitz back¬ 
ground We will now construct our minimal phenomenological model of 

holographic superconductivity by adding to (12.31) the action term 

Sm = j d^+^x^ - \D-^f - m2 |4-|2^ , (2.8) 

where we have introduced a charged scalar field 'k and a C/(l) gauge field following 
[5], and where and — iA^. We will assume that there 

is neglible interaction with the gravitational background and therefore it remains fixed 
and given by the Lifshitz black hole solution (l^-(l^. This lack of back reaction 
means we are effectively working in the probe limit (very large scalar field charge). As 
we will explain below, the scalar field mass will be set so as to get particular dimensions 
for the condensate under study. 

The general equations of motion for these fields are 


DH = m^^-, (2.9) 

= r + \^!\^ A'^, ( 2 . 10 ) 


where 

We propose the following ansatz for the component fields 


^'(m) 


F^(u) 


A = (j){u)dt , 


( 2 . 11 ) 


( 2 . 12 ) 
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where V'(w) is a real function. Under this ansatz the equations of motion (I2.9p - (l2.10p 
become 


V'" + 


f_ d+z-l 

f u 


-2 


1 

U^f 


— 




d — z — 1 , 

u v?f 


V' = 0, 

= 0 . 


(2.13) 

(2.14) 


This system of equations admits the no-hair solution ^(r) = 0. In this case the gauge 
field has solutions 


= p P d-z ^ 

(z^d), 

(2.15) 

(t>{u) = /r-plog^^V 

{z = d), 

(2.16) 


where ^ is a constant. This no-hair solution will correspond to the normal phase of the 
superconductor. The superconducting phase will be given by solutions with ip{u) ^ 0. 
From the equations of motion (I2.13I) - (I2.14I) we see that the asymptotic u —>■ 0 behavior 
of the fields is 


u“- 1I.^+ 


%p{u) Ki 0-—^ -1- 0-f^— + 


(2.17) 


and 


,d—z 


(piu) fV 

^ Pd-z+---^ 

' h 


(2.18) 

(j){u) K. 

/r-plog^^^ +••• , 

iz = d), 

(2.19) 

A± = 

= \[{z + d)±^/{z + d)2 

+ 4771^^ , 

(2.20) 


with 


from where we get the BF-bound on the mass 


( 2 , 21 ) 


According to the AdS/CFT dictionary, the asymptotic coefficient 0+ corresponds to 
an operator of dimension A_|_, while 0_ corresponds to a source in the boundary theory. 
Meanwhile, /i and p correspond to the chemical potential and the charge density of the 
dual field theory, respectively. In order to solve Eqs. (12.13^ - 112.141) we will impose the 
regularity condition at the horizon (j){u = 1) = 0. Also, from Eq. (12.131) we obtain the 
additional condition at u = 1 

V^'(l) = ^^(l)- (2.22) 
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Additionally, in order to simplify the numerical calculations, we will make use of the 
scaling symmetries 

r —>■ ar , < —>■— t, Xi ^—Xi, g^a^g, (2.23) 

a 

As explained in the Introduction, we will focus on the particular case D = 5. This 
means we will take d = 3. When looking for hairy solutions to the equations of motion 
one has two possible boundary conditions at m —>■ 0. One can set 0_ = 0 {standard 
quantization condition), in which case one has the operator 0+ (dimension A_|_) as the 
superconductor order parameter. Conversely, one can set the boundary condition 0+ = 
0 {alternative quantization condition), which results in the operator 0_ (dimension A_) 
being the superconductor order parameter. Having set either one of these boundary 
conditions, one can proceed to solve the equations of motion through the shooting 
method. 

2.2 Different Cases of Condensation 

Going back to the allowed values for the dynamical critical exponent, we see that for 
the D = 5 case we can have 1 < z < 3. Throughout this paper, for both brevity and 
simplicity, we will choose to work with the integer values z = 1,2. This suits perfectly 
our primary objective, stated in the introduction, which is to have a general idea of how 
the dynamical critical exponent z affects our holographic superconductor with respect 
to its behavior in the usual (z = 1) isotropic realization of the gauge/gravity duality. 
As will be seen in the following, the general tendency in the behavior of the physical 
quantities of the system will be very clear when treating these values. 

In order to have a more comprehensive study of the effect of the dynamical critical 
exponent z on our holographic superconductor, we will choose to work in the following 
cases: 

• Case I. We set the value of the scalar field mass as 

m^=-Sz. (2.24) 

In this way, we have 

A_=z, A+=3, (2.25) 

so that the asymptotic behavior of the scalar field at m —>■ 0 is 

'4’{u) ~ Oj,— + O 3 — + • • • . (2.26) 
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Table 1: Value of critical temperature Tc, for different cases. 



z = 1 

z = 2 

Case I 

0.198 

0.087 

Case II 

0.517 

0.351 


In this case, we will set 0^ = 0 for all values of z considered, so that the superconducting 
order parameter of the system will be given by O 3 of dimension 30 

• Case II. We set the scalar mass as 

m^ = -{z + 2). (2.27) 

This choice of mass results in 

A_ = l, A+=z + 2, (2.28) 

so that near u —>■ 0 we have 

U 

tp{u) Ri 0i-h 0^+2 ,2 + • • • . (2.29) 

Th r + 

Here we will choose to set Oz +2 = 0 and the order parameter of the superconductor will 
be given by 0 i of dimension 1 . 

In figures [T][2] we show the value of the condensate Oa as a function of temperature 
for each one of the cases described above. We notice that the near-Tc the condensate 
behaves as 

Oa (1 - T/Te)^/^ , (2.30) 

for all values of z. Therefore, the dynamical critical exponent does not alter the mean- 
field theory behavior of the order parameter. In Table [T] we show the value of the 
critical temperature for our different cases. We notice that the value of the critical 
temperature decreases with z for all cases, and therefore a large dynamical critical 
exponent inhibits superconduction. 


^The same condensate was used in m, but magnetic properties were not studied in that paper. 









Value of the condensate Oa as a function of temperature, for different cases. 



p 


1/3 


P 


2/3 


(a) 2 = 1 


(b) 2 = 2 


Figure 1: Case I. 



pl/3 p2/3 


(C) 2 = 1 


(d) 2 = 2 


Figure 2: Case II. 
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3 Field Fluctuations 


3.1 Gauge Field Fluctuation 


In this section we will add small fluctuations to the component fields of our model. 
As explained before, we will set d = 3. We begin by adding the following gauge field 
fluctuation 

A = (p(u)dt + SAx{t,u,y)dx, (3-1) 


where 

SA,{t,u,y) = e-^^*+^’^yA,{u). 

The corresponding equation of motion for Ax{u) is, to linear order 
d — z — 3 




f 


AZ + 


„22 f2 

'h J 


rlf 


A, =0. 


(3.2) 


(3.3) 


We will consider the case where w and k are much smaller than the scale of the conden¬ 
sate (low-frequency/small-momentum regime), so that quadratic terms in fc, oj in (13. 3 1) 
can be neglected. Demanding regularity at the horizon u = 1, from (1531) we have the 
following conditions 

Aa;(l) = Aa;0, AZ{1) = jy^AxO . (3.4) 

where </'(!) = tpo- Since Eq. (13.31) is linear, we will set A^o = 1 without loss of generality. 
From (13.31) we see that the asymptotic behavior of at u —>■ 0 is 

.d-\-z — 2 


Ax{u) ~ A^*^^ -|- Jx 


d-\-z — 2 


+ 


(3.5) 


According to the AdS/CFT dictionary, aZ^'^ corresponds to a vector potential in the 
dual field theory, while Jx corresponds to its conjugate current [5]. We can relate both 
physical quantities through the London equation 


Jx = -n.Af , 

ms 


(3.6) 


where Ug is the superconducting carrier density number and TOs is the superconductor 
carrier mass. For simplicity we define the quantity 


1 


TOs 


(3.7) 


which can be computed holographically from dSH) and ()3.6I) as fis = —Jx/A^\ In 
figures [3J3] we show the value of hg as function of temperature, for different cases. We 
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Table 2: Value of C'z for different cases. 


a/p(2A-.-l)/3 

z = 1 

z = 2 

Case I 

8.272 

2.068 

Case II 

0.297 

0.032 


find that hg behaves near-Tc as 

^ (1 - T/T,) . (3.8) 

It is found numerically that the ratio of O^/hg as a function of temperature behaves 
almost constantly and has a definite value at T = that varies according to the value 
of 0 within a specific case of condensation. We define this ratio at Tc as 

((i2 

-^ = Cz. (3.9) 

Us 

We show in Table [5] how the constant Cz varies for different cases. In figure [S] we 
show this ratio as a function of temperature, for Case I, z = 1, and Case II, z = 1. 
We can observe that the ratio behaves almost like a constant with respect to T. The 
ratio (1^ will be important in the next section, when we apply the Ginzburg-Landau 
interpretation to our system. 


3.2 Scalar Field Fluctuation 


We now consider a small fluctuation to the scalar field of the form 


with 


+ Si:{u,y)) 


6ij;{u, y) = e*^^? 7 (M). 


The corresponding equation of motion is 




f jd + z-l) 

f u 




I 

v?f 


2 

m — 




ry = 0, 


(3.10) 


(3.11) 


(3.12) 


where we set d = 3. Demanding regularity at the horizon u = 1, from (13.1211 we have 
the following conditions 

'„ 2 ' 


r]{l) = rio, d'(l) = do, 


(3.13) 
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Value of hs as a function of temperature, for different cases. 




(a) a = 1 (b) z = 2 


Figure 3: Case I. 



(c) 2 = 1 {d) z = 2 

Figure 4: Case II. 
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2 
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0.0 
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0.05 0.10 0.15 

T 


0.30 0.35 0.40 0.45 0.50 

T 


P 


1/3 


P 


1/3 


(a) Case 1. z = 1 


(b) Case II. « = 1 


Figure 5: Value of the ratio 0\/fis as a function of temperature, for different cases. 


while at u —>■ 0 we have the asymptotic behavior 

/7(u) (d0_) ^ + (d0+) ^7 + ■ ■ • (3-14) 

'^h 

When solving equation (13.1211 we set the same boundary conditions at u —>■ 0 as for the 
field Ip. Since we will not be concerned about the normalization of rj, we set rjQ = 1. 

Following [22] , we can compute holographically the correlation length of the bound¬ 
ary operator by calculating the wave number k. Indeed, the correlation length is the 
inverse of the pole of the correlation function of the order parameter written in Fourier 
space 

To obtain the wave number k, one must solve Eq. (13.121) as an eigenvalue problem 
consistent with the boundary conditions. This is done near the critical temperature. 
Once having computed k, one obtains the correlation length simply as 


&I = T. (3.16) 

In figures [BjlTI we show the value of fc as a value of temperature for our different cases. 
Also, in figures IHIini we show the value of as a function of temperature, for our cases. 
We find that near the critical temperature. A: ^ (1 — T/Tc)^^^, and equivalently 


^0 


1 

(1-T/T,)^/" ’ 


(3.17) 


for all values of z, which is in agreement with mean field theory. 
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Value of the wave number fc as a function of temperature, for different 

cases. 



{&) z = 1 


(b) z = 2 


Figure 6: Case I. 



(c) 2 = 1 


{d) z = 2 


Figure 7: Case II. 
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Value of the correlation length as a function of temperature, for different 

cases. 




(a) 2 ; = 1 (b) 2 = 2 

Figure 8: Case I. 




(c) 2 = 1 (d) 2 = 2 


Figure 9: Case II. 
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4 Ginzburg-Landau Approach 

At this point we implement a phenomenological Ginzburg-Landau approach to our holo¬ 
graphic superconductor, following m- The main assumption of this approach is that 
the dual field theory can be described near the critical temperature by the effective 
action 

5eff « J dx^ |a + I I'I^glI" + ^ + • • • I . (4.1) 

where the component fields are a scalar field 4 'gl representing the order parameter of the 
theory, and a vector field with ^ = 0,..., 3 and where Di = di — iAi. Also, rus is the 
superconductor carrier mass and a, /3 are phenomenological parameters with tempera¬ 
ture dependence. According to the AdS/CFT dictionary, the vector field components 
Ao and Ax correspond to the chemical potential p, in (12.1811 and to vector potential 
in (13. Sp respectively. According to mean field theory, the order parameter |'I'gl| 
has critical exponent 1/2. In order to match this critical exponent with that of Oa we 
propose the identification 

|4'GL|'=iV,0i, (4.2) 

where Nz is a proportionality constant that depends on z and changes according to 
every model we consider. 

In the remaining of this paper, we adopt the same notation and conventions of |19] . 
In particular, the superconducting carrier mass iris can be absorbed in definitions of the 
other parameters in Ginzburg-Landau theory, and we can therefore safely set nis = 1. 
Going back to (ITtI) . this means in particular that hg = Ug and the numerical equality 
(13.91) can be written as 

(0 2 

^ = G, . (4.3) 

rig 

However, one has according to Ginzburg-Landau theory the following relation between 
the order parameter [TglI and the charge carrier density rig 

|^GL|" = n«. (4.4) 

Then, substituting (1131) and our identification dUl) in (031) we obtain 

A^. = ^ . (4.5) 

In Table HI we show the value of the proportionality constant Nz for various cases. 
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Table 3: Value of Nz for different cases. 


p(2A-z-2)/3^^ 

z = 1 

z = 2 

Case I 

0.121 

0.484 

Case II 

3.367 

30.986 


We can also calculate the penetration length A of the superconductor, defined as 

1 


A = 


\/4t 


(4.6) 


In figures fTHlITT] we show the value of A as a function of temperature, for our different 
cases. As in the case of we have that the behavior of A at T Ri Tc is 

1 


A 


(4.7) 


(1 - T/rc)i/2 ’ 

for all 2 . This result is in agreement with mean field theory. 

In order to have a consistent Ginzburg-Landau description of the dual field theory, we 


must be able to determine by holographic methods the parameters |a| and jS. 


|a|, we can determine it directly from the Ginzburg-Landau theory relatiorjf 


Regarding 


a = 


4 ^ 0 ^ 


(4.9) 


In figures 11^131 we show the value of |a| as a function of temperature, for our various 
cases. We see that the near-T, behavior of |a| is 


|a| ^ ao (1 - T/T,) , 


(4.10) 


which is in agreement with usual Ginzburg-Landau theory, for all 0 . However, we find 
numerically that the value of the coefficient ag decreases as the value of z raises. 

The remaining phenomenological parameter {3 can be computed through the Ginzburg- 


Landau theory relation 




(4.11) 


^The actual Ginzburg-Landau theory relation is 


i«i = 




( 4 . 8 ) 


where ^ is the Ginzburg-Landau coherence length. The coherence length is in turn related to the 
correlation length ^ = 2 ^q. See m- 
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Value of the penetration length A as a function of temperature, for 

different cases. 

10 
8 
6 
4 
2 
0 

0.05 0.10 0.15 0.20 

T 

{&) z = 1 

Figure 10: Case I. 

1.2 

1.0 

0.8 

0.6 

0.4 

0.2 

0 . 0 ' 

0 . 



{d) z = 2 


2 - 



p'/^A 




(h) z = 2 


Figure 11: Case II. 
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Value of the parameter a as a function of temperature, for different cases. 




(a) 2 = 1 (b) 2 = 2 

Figure 12: Case I. 




(C) 2 = 1 


(d) 2 = 2 


Figure 13: Case II. 
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Table 4: Value of the Ginzburg-Landau parameter k, for different cases. 


n 

z = 1 

z = 2 

Case I 

0.527 

0.467 

Case II 

0.070 

0.002 


where |4'oo| is the value of the condensate deep inside the superconductor, where external 
fields and gradients are negligible. Since we are in the limit of small field perturbations, 
we indeed find ourselves in that approximation. Substituting (03) and (14.9|) in (I4.11|) 
we obtain the following expression 


/3 = 


1 1 


(4.12) 


In figures [TUllSI we show the behavior of /3 as a function of temperature, for our different 
condensation cases. We observe that, near-Tc, (3 behaves in agreement with Ginzburg- 
Landau theory, having a definite value at T = T^- We also observe that this value 
decreases as the value of z raises. 

Finally, having calculated the characteristic lengths of the system and A, we can 
compute the Ginzburg-Landau parameter k, defined as 


A 



(4.13) 


where ^ is the Ginzburg-Landau coherence length, which is related to our correlation 
length ^0 as = 2^g. (See [IH])- In figures [TH1IT7I we show how the Ginzburg-Landau 
parameter n behaves as a function of temperature, for our different cases. We notice 
that all plots have a definite value at T = Tc. We will take this to be the value of k of 
our holographic superconductor for each case considered. The value of k for different 
cases is shown in Table 01 We note that all values of k are lower than l/-\/2 ^ 0.707 
for all cases of z considered, which means that our system behaves always as a Type I 
superconductor. Also, we notice that the value of k is always lower for z = 2, which 
means that in holographic superconductors with higher dynamical critical exponent, 
vortex formation is more strongly unfavored energetically and have a stronger Type I 
behavior. 
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Value of the parameter /3 as a function of temperature, for different cases. 




(a) 2 = 1 


(b) 2 = 2 


Figure 14: Case I. 




(C) 2 = 1 


(d) 2 = 2 


Figure 15: Case II. 
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Value of the Ginzburg-Landau parameter k as a function of temperature, 

for different cases. 




(a) 2; = 1 (b) 2 = 2 


Figure 16: Case I. 




(c) 2 = 1 (d) 2 = 2 


Figure 17: Case II. 


22 







5 Constant Magnetic Field 


We will now study the effect of an external magnetic field to the superconducting phase 
of our models. As done before, we begin with in the general dimensional case and then 
focus on Z? = 5. We follow the procedure developed by [20] and will proceed in a 
perturbative fashion by proposing a series expansion for the component fields 

^{x,u) = -\ - (5.1) 

A^{x,u) = A^°\x,u)+eAA\x,u)-{ - (5.2) 


where x = {x, y), and the expansion parameter is given by 


Ba-B 


e < 1, 


(5.3) 


were Be is the value of the magnetic field that breaks the superconducting phase (critical 
magnetic field). Since this expansion is done near the value B = B^, this means that we 
find ourselves near the point where the condensate vanishes. We substitute expansions 
(I5.1[) - (I5.2I) in the general equations on motion (I2.9I) - (I2.10I) . The zero order equation for 


the gauge field is 





Va 

(v^doi)=». 

(5.4) 

and has solutions 





A^f\u) 

^d-z 

(5.5) 


Af\x) 

= BeX, 

(5.6) 


and the rest of spatial components equal to zero: = 0, i ^ y. Since the solution for 

is equal to solution (12.151) . we set the notation = (j), for simplicity. Meanwhile, 
the general scalar field equation is 


d-\-l — z 


where I,J = x, y. Eq. dSjJ is clearly separable. We follow the standard treatment and 
propose 

{x, u) = {x, u ), (5.8) 


so on the right hand side of (EB we have 


= (92 + _ iBe x)^) = e*™ -{p-Be x)"^^ , (5.9) 
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and we get the following equation 


- (5r - ;p) (-9- + (f - »)’) ■ 

(5,10) 

Now, we make the separation 


and define the variable 


‘Pn\x,u) = Pn{u)'y^'> {X) , 


X = J ’ 


so that the operator on the right hand side of (I5.10|) becomes 


-d^ + (p-B,xf =(2B,) 


-d\ + 4x2 


iv) 

and acting on we have the eigenvalue equation 

(-aj + , 

that has as a solution the eigenfunctions 


(5.11) 


(5.12) 


(5.13) 


(5.14) 


(5.15) 


with eigenvalues 


A„ = 2n + 1, n = 0,1... 


(5.16) 


We choose the n = 0 mode, which corresponds to the most stable solution Eiiiniiis]- 
As described originally in [20], the more general solution to the scalar field is given 
by linear superposition of the solution obtained above, with different values of p. (We 
adopt the authors notation in the following). Going back to (15.81) . (15.lip and (I5.15p . we 
write our solution explicitly as 


00 

(u,f) = po(w) X! Cie"^‘yjo{x;pi) , 

I —— 00 


where 


and where we define 


7 o {x;pi) = exp 



Cl = exp 


-I- 


2'n:y/Bl 

Pi = - 1 , 

Oi 


(5.17) 


(5.18) 


(5.19) 


24 
















and oi, 02 are real parameters. Solution (15.171) can be rewritten as 


{u,x) = ^po{u)e 2 d 3 (t;,r) 
where i? 3 (r;,r) is the elliptic theta function, defined as 

OO 


e e 


and where the variables v and r are defined as 


1 , 

T = ^ (2i7r — 02 ) . 


(5.20) 


(5.21) 


(5.22) 


Q-i U/^ 

Owing to the elliptical theta function '^ 3 , the scalar field solution has the 
following pseudo-periodicity in the x — y plane 


(u,x,y) = a;, 2 /-I-oi). 


(5.23) 


( u, X -f 


27r 


-,y + 


02 




= exp 




Ol 


02 


2(21 


iu,x,y) . 


In addition to this, the tDs function has zeros located periodically at 


1 


l/= r72-|-- ui-|- n-|--W72, 


1 


where the fi) vectors are given by 
^ ai d 

Vi = 


27r d 


y/Bf dy ’ v^oi dx ' y/Bfa 2 dy 

Thus, the solution has a lattice profile in the [x—y) plane, spanned by the vectors u). 
We note that, in our given approximation, we will get a 2-dimensional plane, orthogonal 
to the remaining (cZ—2)-dimensional boundary space, where the vortices live. We should 
note that the presence of the vortex solutions given above does not contradict the fact 
that our system was found in the previous section to be Type Indeed, the computation 
of K presented above comes from an energetic analysis, conducted directly from the dual 
system’s Ginzburg-Landau action. (See [IH])- This shows that, according to Ginzburg- 
Landau theory, the formation of the above vortex solutions costs more energy to the 
system than the energy needed for the system staying in a superconducting state. (See, 


02 d 


(5.24) 


(5.25) 


(5.26) 


^ A dynamical approach to vortex solutions in Z) = 4 can be found in m, where it was concluded 
that, for some values of the system’s parameters, the dual superconducting system was Type I 
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Value of the critical magnetic field Be as a function of temperature, for 

different cases. 



Be 



0.03 0.04 0.05 0.06 0.07 0.08 

T 


„ 2/3 


{&) z = 1 


(h) z = 2 


Figure 18: Case I. 




(c) z = l 


{d) z = 2 


Figure 19: Case II. 
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Returning to the scalar field equation (15.101) and substituting the results given above, 
we obtain the following equation for the radial function p 


d-\-l — z 


du 


f 


iZ-\-d— 1 


dup{u) I - 


Br 


j2z J.2ZJ j,2^y2z-2 


p{u) = 0 , 


which can be written as 

7' d + z-i 


p 


f 


1 

u^f 


— 




This equation of course has the same behavior at u —>■ 0 as (12.171) 

pr^C-U^- +C+u^+ , 


(5.27) 


(5.28) 


(5.29) 


with A± given by (I2.20|) . We set the same boundary conditions at u —>■ 0 as for the field 
Ip in (I2.12p . By applying the shooting method to Eq. (I5.27P we find the value of the 
critical magnetic field that breaks the superconducting phase of the system. In figures 
mm we show the value of the critical magnetic field Be as a function of temperature, 
for our different cases. We see that near-Tc the critical magnetic field Be behaves as 


Be^{l- T/Te) , 


(5.30) 


which is in agreement with mean field theory, for all values of z. We also note by 
comparing Eqs. (15.281) and (I3.12p . that the procedure to obtain the near-Tc values of 
the square of the wave number k and the critical field Be is the same. This in turn 
confirms the relation between the correlation length and the critical magnetic field put 
forward in m 

Be-^, {T-Te). (5.31) 

so 


6 Conclusions 

In this paper we have constructed a, D = 5 minimal model of holographic superconduc¬ 
tivity in the probe limit, with a Lifshitz black hole background. Within this framework, 
we have studied different cases of condensation, varying within each of them the dy¬ 
namical critical exponent in order to gain insight on how the system is affected by z 
with respect to its usual isotropic behavior. We have added small scalar and gauge field 
fluctuations to the original component fields. These fluctuations allow us to compute 
holographically the penetration and coherence length of the superconducting system. 
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We saw that both characteristic lengths have the standard near-Tc functional depen¬ 
dency on temperature for all condensate cases and all values of z. However, the dynam¬ 
ical critical exponent z does affect the value of the characteristic lengths, as it becomes 
evident in the change of the value of their ratio as given by the Ginzburg-Landau pa¬ 
rameter K. We also saw that it is possible to construct a consistent Ginzburg-Landau 
phenomenological interpretation of the dual theory with Lifshitz scaling. We computed 
through holographic techniques the Ginzburg-Landau Lagrangian parameters a, /3 and, 
as with the characteristic lengths, concluded that they have the standard near-Tc func¬ 
tional dependency on temperature for all condensate cases and all values of z. However, 
the presence of z does have a non-trivial effect on this phenomenological parameters, 
diminishing the value of their numerical coefficients as z raises. 

We have also computed with holographic techniques the Ginzburg-Landau param¬ 
eter K of the system. For all case of condensation and all values of z, we saw that 
K < This means that for all cases the dual system will behave as a Type I super¬ 

conductor. Moreover, we also observed that, for each case of condensation considered, 
the value of n, became lower for higher values of z. This means that in systems with 
higher anisotropy, vortex formation is more strongly unfavored energetically and exhibit 
a stronger Type I behavior. 

Finally, we computed the critical magnetic field needed to break the supercon¬ 
ducting phase of the system, following the perturbative procedure first developed in 
m- We observed that the critical field near-Tc functional dependence on temperature 
is the one predicted by Ginzburg-Landau theory. However, we also note that the value 
of the critical magnetic field is smaller for higher values of z. Additionally, within this 
perturbative approach, we have confirmed holographically the conjecture posed in HZ] 
that the critical magnetic field is inversely proportional to the square of the correlation 
length, in accordance to Ginzburg-Landau theory. 

All of the above results were obtained from a minimal model of superconductivity 
following [Sj. It would be interesting to see how these results would be affected by 
the choice of other models, such as, for instance, d-wave holographic superconductors 
|25) . models with higher corrections to the scalar field potential such as the ones that 
appear in top-down approaches [ZSl HZl HH] or less conventional models such as ones 
with Chern-Simons terms, higher-derivative couplings or in the context of New Massive 
Gravity HU 1301 El- 
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